Magnetic field induced chiral particle-hole condensates 
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We demonstrate that a chiral particle-hole condensate is always induced by a number-conserving 
ground state of non-zero angular momentum in the presence of a magnetic field. The magnetic 
interaction originates from the coupling with the intrinsic orbital moment of the chiral state when 
the field is applied perpendicularly to the plane. According to our numerical results the induction 
mechanism is practically temperature independent providing robustness to these states up to high 
temperatures. This opens the door for manipulating the anomalous Hall response resulting from this 
intricate class of states for technological applications while it also suggests that chiral particle-hole 
condensates may be hidden in various complex materials. 
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Among the numerous states that emerge in a strongly 
correlated system in the particle-hole channel, the chi- 
ral states enjoy special attention. These states apart 
from the related symmetry breaking they effect due to 
the corresponding pair condensation, they additionally 
violate parity and time-reversal because of their special 
momentum structure (See e.g. As a consequence, 

the condensate carries an orbital moment that can be 
viewed as the effect of a non-zero Berry curvature in 
fc-space 0, H, |3l- The arising orbital moment @ pro- 
vides the system with a magnetic field coupling leading 
to anomalous Hall transport and enhanced diamagnetic 
response. As a matter of fact, such states have ideal 
properties for technological applications. In addition the 
well established belief that unconventional particle-hole 
condensates are hidden in the majority of several impor- 
tant materials render these states as strong candidates 
for a pleiad of unidentified phases. Consequently, it is of 
high priority to be in position to generate these states 
and manipulate them, at our own disposal. 

One of the traditional ways to controUably engineer 
novel states of matter is via applying an external field. 
Typical examples of field induced states in the area 
of particle-hole condensates constitute the field induced 
charge and spin density waves. There are several kinds 
of field induced density waves such as those occurring 
in organic quasi-lD conductors Q or even confined spin 
density waves [71 both generated due to the orbital cou- 
pling with the magnetic field. Nevertheless, one may also 
obtain induced density waves driven by the Zeeman cou- 
pling i. 

In this letter we perform a numerical study of the 
magnetic field induced planar chiral particle-hole con- 
densates. A similar situation was first addressed by R. 
B. Laughlin in the context of chiral d-wave supercon- 
ductors to explain the magnetic field induced transition 
observed in Bi2Sr2CaCu208, and later extended to the 
particle-hole channel by J.-X. Zhu et al [l^]. Since the 
whole class of these chiral states are characterized by a 
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universal behaviour we shall concentrate on a particular 
state belonging to this class, the chiral d-density wave 
state, that recently attracted a lot of attention due to 
its prominent role in explaining the pseudogap regime of 
the cuprates 11, 12, 13,, [ij, O [l^. By taking into ac- 
count the effect of an external perpendicular magnetic 
field through its coupling with the intrinsic orbital mag- 
netic moment of this state, we demonstrate through a de- 
tailed numerical analysis that a chiral state is necessarily 
induced or it is strongly enhanced if it already exists. In 
addition, we extract the magnetic field dependence of the 
two density wave order parameters, which is different in 
the above two cases. However, the interplay of the un- 
derlying interaction and the field strength can lead to a 
crossover in the magnetic field dependence in the latter 
case. Furthermore, we observe that the chiral d-density 
wave is robust against the increase of temperature in the 
presence of the external magnetic field, a direct conse- 
quence of the field driven enhancement. Our results are 
striking demonstrating that in many materials in which 
possibly only some components of a chiral particle-hole 
condensate develop, there will he unavoidably an induc- 
tion of the rest and a concomitant transition to the com- 
plete chiral state in the presence of an external field, giv- 
ing rise to the aforementioned intriguing response. 

As we have already mentioned, the representative chi- 
ral particle-hole condensate that we shall consider is 
the chiral d-density wave state that constitutes a sin- 
glet unconventional density wave with a planar momen- 
tum structure, characterized by the commensurate wave- 
vector Q — [t: / a,-K / a). It is composed by a real dxy 
charge density wave violating parity and an imaginary 
dx2-y2 orbital anti-ferromagnetic state giving rise to local 
charge currents and zero charge density, violating time- 
reversal. This state has been shown to exhibit unconven- 
tional Hall transport and anomalous magnetic response 
which is common to any other chiral particle-hole con- 
densate with zero [l7j or finite momentum taking place 
either in the singlet or triplet channel of some kind of a 
spin degree of freedom, that is characterized by the same 
Berry curvature. 

First of all, the chiral d-density wave is known to giv e 
rise to the Spontaneous Quantum Hall effect [l^, U^ . 
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which concerns the generation of a quantized Hall volt- 
age via the sole application of an electric field. Quite 
similarly, a thermoelectric Hall effect can be reproduced 
by the application of a finite temperature gradient T^]. 
Moreover, this state supports topological spin transport 
characterized by dissipationless spin currents in the pres- 
ence of a Zeeman field gradient p^ . An even more 
striking behaviour dominates the magnetic response of 
this system. The existence of the intrinsic orbital mag- 
netic moment leads to perfect diamagnetism and conse- 
quently to the Topological Meissner effect demonstrated 
in Ref. [l^ . In the half- filled case the topological Meissner 
effect is identical to the usual superconducting diamag- 
netism, and motivated by this feature the chiral d-density 
wave state has been proposed to be hidden in the pseu- 
dogap phase of the cuprates. Furthermore, the observed 
Polar Kerr effect in YBa2Cu306+2; [19], has been con- 
sidered as a sharp signature of the chiral d-density wave 
state in underdoped cuprates [ll|. 

To model the chiral d-density wave state we consider 
the following Hamiltonian 
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where we have introduced the single band energy disper- 
sion of the free Bloch electrons e(fe) = —2t[cos{kxa) + 
cos{kya)] arising from the nearest neighbours hopping 
term and a 4-fermion interaction driven by an effec- 
tive separable potential V(fc, fc') = Vifi{k)fi{k') + 
V2f2{k)f2{k'), projecting only onto the form factors 
/i = sm{kxa) sin(fcj,a) and /2 — cos{kxa) — cos{kya) of 
the dxy and dx^_y2 momentum space orbitals. The oper- 
ators Cfc/Cfc annihilate/create an electron of momentum k 
in the reduced Brillouin zone (R.B.Z.), v is the volume of 
the system and a the square lattice constant. For conve- 
nience we have omitted the spin indices, as we are dealing 
with a spin singlet state. Within a mean field decoupling 
we obtain the order parameter of the chiral d-density 
wave state A(fc) = Ai sm{kxa) sa\{kya) — il^2[cos{kxa) — 
cos(fcya)], which satisfies the self-consistence equation 

^(^) = -iE^(fc'^')(4'Cfc-+g), (2) 
k' 

with <> denoting thermal and quantum-mechanical av- 
erage. For convenience, we adopt a more compact nota- 
tion by considering the enlarged spinor ^E"^ = {c\. c^k^o) 
and employing the t Pauli matrices. The mean field 
Hamiltonian now becomes H = T.k'^l'^i^)^ k = 
^fc 9{^) ■ ^ k where we have introduced the isovec- 
tor g{k) defined as g{k) = {TZeA{k),—TrnA{k),e{k)). 
The one-particle Hamitonian for each fc-mode is a 2 x 2 
matrix 



n{k) = 
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33 (fc) -*.92(fc) 

+i52(fc), 753(fc) 
To incorporate the interaction with the magnetic field 

we shall consider only the orbital coupling and neglect 



the Zeeman term as it is usually negligible in the case 
we are considering. To calculate the orbital moment one 
has to determine the Berry phase emerging in this chi- 
ral state when k changes adiabatically along a closed 
loop. For instance, the variation of k can be enforced 
by the minimal coupling k ^ k e£t if we apply 
a constant electric field £. In this case, the Hamil- 
tonian becomes parametric H{k) — > 7i{k,t). Within 
the adiabatic approximation, the emerging Berry phase 
can be determined using the instantaneous (snapshot) 
eigenstates of the parametric Hamiltonian |<E>^(A;,t), sat- 
isfying n{k,t)\^^{k,t)) = E^{k,t)\^^{k,t)). In this 
equation, is introduced only as a parameter. This 
means that these snapshot eigenstates are not really 
time-dependent but only parameter dependent, which 
in our case coincides with t. Our two band system, 
is characterized by the snapshot eigenstates |$±) and 
the corresponding eigenenergies E±{k,t) — ±\g{k,t)\. 
By defining gi{k) = E{k) sin 9 (k) cos (p{k), g2{k) — 
E{k) sin 6 (k) sm If (k) and g^ik) = E{k) cos 9{k), we ob- 
tain a convenient expression for the snapshot eigenstates 
of the system 
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with E{k) — \g{k)\ and denoting matrix transposition. 
The orbital magnetic moment mj^{k) is determined by 
the relation m,(fe) = ^(V^a>,(fc) | x [n{k) ~ E,{k)] \ 

V^<i>^(fc)) = ^E{k)Cl^{k) where we have introduced 
the Berry curvature [15] 
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We notice the Berry curvature and the orbital moment 
lie along the z-axis, as a direct consequence of the planar 
character of our system. The presence of a perpendicu- 
lar magnetic field B enters the band dispersions of the 
system in the following way E^{k) = Ej^{k) — m^ik) ■ B. 

Having obtained the energy dispersions of the system 
in presence of the magnetic field we may now extract the 
self-consistence equations of the dxy and dx^-y'^ order 
parameters Ai, A2. The free energy functional is defined 
as 



_ A2 A2 1 V- / 
= — — - — > Ln(l + e- 
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with P — \/kgT. The first two terms in the free en- 
ergy originate from the mean field decoupling and corre- 
spond to the elastic energy waisted for building up the 
two density wave gaps. Minimization of this functional 
with respect to the order parameter doublet leads to the 
following system of coupled self-consistence equations 
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z = l,2, (7) 



that will be used for determining the chiral d-density 
wave gaps numerically. In the above, we have introduced 
s{k) = sin^(fca;a) cos^(fcya) + sin^(fcj^a), Ai^2 = A2,i, 
np{E) the Fermi-Dirac distribution and multiplied by 
a factor of two in order to take into account the electron 
spin. The first term corresponds to the equation that we 
obtain in a zero magnetic field with the only difference 
that the energy bands are shifted by the orbital cou- 
pling ■ 3. On the other hand, the second term is 
attributed entirely to the magnetic field interaction with 
the chiral d-density wave state. As a matter of fact, it is 
the essential ingredient for a field induced chiral d-density 
wave. To understand how a magnetic field stabilizes this 
chiral state we consider the case of Ai = and A2 ^ 0. 
In this case the initial state is an orbital anti-ferromagnet. 
We may calculate the induced Ai component by setting 
Ai = on the right hand side of Eq.([7]). It is straightfor- 
ward to obtain that = SS^B^V^ ^2^(6,) with 
^{B^) corresponding to a suitable sum over R.B.Z. of the 
remaining terms. Apart from the additional weak depen- 
dence on the field in T, this formula also agrees with the 
one derived in Ref.[10] using a different method. Accord- 
ing to our result a chiral d-density wave state is always 
generated even for an arbitrary small magnetic field. 

We now proceed in solving numerically the system of 
self-consistence equations. For all our numerical simula- 
tions we set t — 250meV, a — 5 A and V2 — ISOmeV. The 
last condition aims to establish a dx2-y2 density wave of 
approximately a A2 ~ 53meV gap for all the possible 
temperatures and magnetic fields considered here. The 
latter consideration helps us to focus solely on the be- 
haviour of the dxy component. For the calculations we 
have set up an 128 x 128 grid in the right upper quad- 
rant of the Brillouin zone. A large number of points are 
needed in order to stabilize a value for the d^y order pa- 
rameter. 

First of all, we verify numerically the linear depen- 
dence on the magnetic field of the induced dxy or- 
der parameter. For the illustration we consider Vi — 
125meV, a value that in zero magnetic field would 
provide a practically zero dxy gap- In the inset of 
Figd] we do observe the expected linear scaling while 
we also notice the extremely weak temperature depen- 
dence. The latter is attributed to the manner in which 
the Fermi-Dirac occupation numbers enter the second 
part of the self-consistence equations. Specifically, the 
occupation numbers of the two bands add up giving 
Up. [-E{k) - m^{k)B^] +np [E{k) - m^{k)B^] ~ 1 due 
to the negligibleness of the magnetic coupling compared 
to E{k). 
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FIG. 1: (Color online) Magnetic field induced chiral d-density 
wave. In the inset, we present the magnetic field dependence 
of the induced dxy component when only a dx2_y2 already 
exists, for two diflterent temperatures. In the main panel we 
demonstrate the strengthening due to the field, of an initially 
subdominant d^y component of a preexisting chiral d-density 
wave state. The two curves presented for the dxy part differ 
only on the zero field gap value arising from the usual BCS 
temperature dependence. As we may observe, the induced 
part of the chiral state is always temperature insensitive, while 
the field dependence changes from linear to square root. 



More intriguing results emerge in the case of an ini- 
tially existing chiral d-density wave state. We con- 
sider the combination of potentials Vi — 550meV and 
V2 = ISOmeV, generating the gaps Ai ~ l.SmeV and 
A2 ~ 53meV in the absence of the field. When the or- 
bital magnetic coupling is switched on, the dx^-yi gap 
remains mostly unaffected while the dxy component is 
significantly enhanced. Specifically, the magnitude of the 
gap becomes about 15 times bigger for B^=^QT. More- 
over, we observe that the magnetic field dependence of 
the latter gap turns out to be square root contrary to the 
linear dependence obtained earlier. Of course, this dif- 
ference arises from the first term of the self-consistence 
equation which dominates the zero field limit and is fully 
active here, compared to the previous situation. 

As far as the temperature dependence is concerned, 
we notice in the main panel of Fig[l] that for two to- 
tally different temperature regimes we only have a shift 
of the two curves, equal to the difference of the dxy gap 
value obtained in these temperatures in the absence of 
the field. Apparently, the induced part of the dxy gap 
is once again temperature independent as it originates 
solely from the magnetic field coupling term. We also 
obtain the temperature evolution for different magnetic 
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FIG. 2: (Color online) Magnetic field dependence of the in- 
duced part of the d^y component in a preexisting chiral d- 
density wave state for several interaction potentials. For small 
values of the potential we obtain the expected square root 
behaviour, while after a critical value a crossover to a linear 
dependence is triggered. The latter, signals the dominance of 
the magnetic field coupling over the zero field 4-fermion d^y- 
channel interaction strength. (For the calculations we have 
used V2 = 150meV.) 

field values. For example a Ai = l.SmeV gap initially 
disappearing at about 30K is now robust over a temper- 
ature range of more than 120K. This is natural if we take 
into account that the magnetic field has strengthened the 
zero temperature gap of the order parameter that would 
now collapse to a higher temperature following the usual 
BCS behaviour. 

Finally, we examine the influence of the d^y interac- 
tion potential, on the corresponding gap magnetic field 
dependence. Fig[2] shows the field dependence of the dxy 
gaps for several potentials after subtracting the zero field 
contribution. Focusing on the induced part of the gap 
we conclude that the increase of the interaction strength 
softens the square root field dependence turning it into a 



linear one. This crossover behaviour is indicative of the 
enhancement of the magnetic field coupling that favours 
a linear trend. As a matter of fact, this features serve as 
a potential diagnostic method for the interaction energy 
scale. 

In conclusion, we have studied numerically the occur- 
rence of magnetic field induced chiral particle-hole con- 
densates due to the coupling of their intrinsic orbital mo- 
ment to a perpendicular magnetic field. According to 
our results, if an incomplete chiral state develops in zero 
field, such as a dx2_y2 density wave, then a d^y is gen- 
erated spontaneously when the magnetic interaction is 
triggered. The magnetic field dependence of the dxy gap 
is linear. On the other hand, if a small dxy compared 
to the dx2-y2 already exists, then it will have a square 
root dependence on the magnetic field that can become 
linear if the interaction potential in the dxy channel ex- 
ceeds a critical value. In both cases we obtain a negligible 
dependence of the induced order parameters on tempera- 
ture. Consequently, chiral particle-hole condensates can 
survive up to high temperatures as long as the magnetic 
coupling persists. 

The unavoidable transition to a chiral state has strong 
impact on materials that are proposed to host orbital 
anti-ferromagnetic states. For instance, the well estab- 
lished d-density wave scenario [1^ in the pseudogap 
regime of the cuprates provides a unique occasion for the 
realization of a chiral state, the chiral d-density wave. In 
this case, its presence is ensured when an external per- 
pendicular field is applied but it could also be present in 
zero field if the necessary orbital interaction originates 
from intrinsic sources such as magnetic impurities. 
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